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Nonlinear Motion in Electromagnetic Fields via Multiscale Expansions 

The consideration of dynamics of relativistic beams/particles is based on variational approach to rational (in dy- 
namical variables) approximation for equations of motions. It allows to control contribution from each scale of 
underlying multiscales and represent solutions via exact nonlinear eigenmodes expansions. Our approach is based on 
methods provided possibility to work with well-localized bases in phase space and good convergence properties of the 
corresponding expansions. 

We consider as the main example the particle motion in storage rings in standard approach based on consideration in 
[1]. Starting from Hamiltonian, which describes classical dynamics in storage rings Ti(f, P, t) — c{n 2 + m^c 2 } 1 I 2 + e<j) 
and using Serret-Frenet parametrization, we have after standard manipulations with truncation of power series 
expansion of square root the rational approximation (up to n-poles) for Hamiltonian of orbital motion in machine 
coordinates. So, our problems may be formulated as the systems of ordinary differential equations 

dx; 

Qi( x )~3T = Pi{ x it)i X = (xi,...,x n ), i = l,...,n, ma,xdegPi=p 7 maxdeg Qi = q (1) 
at ii 

with initial (or boundary) conditions Xi(0), Xi(T) and Pi,Qi are not more than polynomial functions of dynamical 
variables Xj and have arbitrary dependence on time. Of course, we consider such Qi(x) which do not lead to the 
singular problem with Qi(x), when t = or t = T, i.e. Qi(x(0)), Qi(x(T)) ^ oo, 0. We'll consider these equations as 
the following operator equation. Let L be an arbitrary nonlinear (rational) matrix differential operator of the first 
order with matrix dimension d corresponding to the system of equations (1), which acts on some set of functions 

* = V(t) = (V (£),..., * d (f)Y teflcR from L 2 {Vt): 



= L(R,t)^[t) = 0, R = R(t,d/dt,V). (2) 

Let us consider now the N mode approximation for solution as the following ansatz (in the same way we may consider 
different ansatzes): 

N 

= $>?ty r (*) (3) 

r=l 

We shall determine the coefficients of expansion from the following variational conditions (different related variational 
approaches are considered in [l]-[4]): 

L% = J (LH> N )Mt)dt = (4) 

We have exactly dN algebraical equations for dN unknowns a r . So, variational approach reduced the initial problem 
(1) to the problem of solution of functional equations at the first stage and some algebraical problems at the second 
stage. Here 4>k(t) are useful basis functions of some functional space (L 2 , L p , Sobolev, etc) corresponding to concrete 
problem. As result we have the following reduced algebraical system of equations (RSAE) on the set of unknown 
coefficients af of expansions (3): 

LiQi^af , ai ) = M{P ih a?,pj), (5) 

where operators L and M are algebraization of RHS and LHS of initial problem (|l|) . Qij are the coefficients (with 
possible time dependence) of LHS of initial system of differential equations ([!]) and as consequence are coefficients 
of RSAE. Pij are the coefficients (with possible time dependence) of RHS of initial system of differential equations 
(Q) and as consequence are the coefficients of RSAE. I = (i l5 i q +2), J = (ji, •••jj'p+i) are multiindexes, by which 
are labelled ctj and /3j, the other coefficients of RSAE ([|): 

Pj = = / II 011 = K-^+J = X] /vv-^-^v^ ( 6 ) 

l<ifc<P+l ill- 
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where p (q) is the degree of polynomial operator P(Q) (0), ie — (1, q + 2), ipi a — d^/dt. According to [l]-[4] we 
may extend our approach to the case when we have additional constraints on the set of our dynamical variables \& 
or x. In this case by using the method of Lagrangian multipliers we again may apply the same approach but for the 
extended set of variables. As result we receive the expanded system of algebraical equations analogous to the system 
(5). Then, after reduction we again can extract from its solution the coefficients of expansion (3). Now, when we 
solve RSAE (||) and determine unknown coefficients from formal expansion (3) we therefore obtain the solution of 
our initial problem. It should be noted if we consider only truncated expansion (3) with N terms then we have from 
(||) the system of N x d algebraical equations with degree £ = max{p, q} and the degree of this algebraical system 
coincides with the degree of initial differential system. So, we have the solution of the initial nonlinear (rational) 
problem in the form 



JV 



x(t)=x(0) + Y,akMt), (7) 

fe=i 

where coefficients are roots of the corresponding reduced algebraical (polynomial) problem RSAE (|J). Conse- 
quently, we have a parametrization of solution of initial problem by solution of reduced algebraical problem (^|). 
The problem of computations of coefficients a] , f3j (^|) of reduced algebraical system may be explicitly solved 
in wavelet approach. The obtained solutions are given in the form (|t]), where ipk{t) are wavelet basis functions. 
In our case ipk{t) are obtained via multiresolution expansions and represented by compactly supported wavelets. 
Because affine group of translation and dilations is inside the approach, this method resembles the action of a mi- 
croscope. We have contribution to final result from each scale of resolution from the whole infinite scale of spaces: 
...V-2 C V—i C Vo C V\ C V% C where the closed subspace Vj(j £ Z) corresponds to level j of resolution, 
or to scale j. This multiresolution functional space decomposition corresponds to exact nonlinear eigenmode de- 
compositions (7). It should be noted that such representations give the best possible localization properties in the 
corresponding (phase)space/time coordinates. In contrast with different approaches formulae (7) do not use pertur- 
bation technique or linearization procedures and represent dynamics via generalized nonlinear localized eigcnmodes 
expansion. So, by using wavelet bases with their good (phase)space/time localization properties we can construct 
high-localized (coherent) structures in nonlinear systems with collective/complex behaviour. As a result our N mode 
construction (7) gives the following representation for solution of equations (1): 

x(t) = xt w (t) + J2^M, "«~2* (8) 

i>N 

where x r (t) may be represented by some family of (nonlinear) eigenmodes and gives as a result the multiresolu- 
tion/multiscale representation in the high-localized wavelet bases. 
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